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This paper presents a sizing optimization technique for Hybrid Stand-Alone 
Photovoltaic (HSAPV). In this research, three optimization techniques have 
been developed, namely Dolphin Echolocation Algorithm (DEA), East 
Evolutionary Programming (PEP), and Classical Evolutionary Programming 
(CEP). These techniques have been incorporated into the sizing process to 
maximize the technical performance of the SAPV system. The components 
of PV modules, charge controllers, inverters, and batteries are used to 
determine the optimum value. These components are used as the control 
parameters to maximize the expected performance ratio (PR) of the SAPV 
system. The Iterative Sizing Algorithm (ISA) is the benchmarking technique 
to conduct the optimization technique achieving maximum PR value and 
minimal computation time. Results obtained from the research show that DE 
overcomes FEP and CEP. In addition, the optimization techniques also 
demonstrated comparatively fast with respect to ISA as the benchmark 
technique. 
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I. INTRODUCTION 

Solar energy is one of the promising types of renewable energy resources in Malaysia with a 
potential total installed solar electricity capacity of more than 6500 MW [1]. The solar electicity is generated 
using either the Grid-Connected Photovoltaic (GCPV) systems or Stand-Alone Photovoltaic (SAPV) 
systems. Mostly, GCPV systems which are widely installed in urban areas where the supply from utility grid 
is highly available. An SAPV system typically consists of PV arrays, solar batteries as well as power 
conditioning units such as charge contollers and inverters. Besides that, the system could also be hybridized 
with other energy resources such as wind, hydro and diesel generator. Moreover, the type, capacity, size and 
configurations of these system components are commonly selected through a systematic process known as 
sizing [1], [2]. Eor PV-battery-diesel generator normally called as a HSAPV systems. 

HSAPV system usually consists of PV modules, such as power conditioning unit charge controller 
and inverter, a battery bank for storage charges and diesel generator (DG). The size of the HSAPV system 
appears to be more complex and difficult than the size of the photovoltaic system grid-connected since more 
components should be considered in the HSAPV system. In [3], [4], the sizing parameters of the HSAPV 
system is more complicated since the system cannot receive immediate back-up from the utility grid. The 
size of the HSAPV system usually requires system designers to initially choose PV modules, charge 
controller, battery bank, inverter and diesel generator before trying to match the electrical characteristics 
among these components. 
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However, if there is more than one HSAPV model, the sizing process is repeated to achieve the 
most suitable parameters for the design. The conventional sizing process has experienced computational 
burden as reported in [5], [6]. Several studies have been conducted to accelerate the overall sizing process in 
HSAPV design. Most approaches made use of various types of Computational Intelligence (Cl) techniques. 
One of the popular optimization techniques that have been widely employed is Artificial Intelligence (AI) 
[7]. AI technique has dealt with difficult problems in some areas; which in turns alleviated the setback 
experienced in the existing conventional techniques. 

There are three important components in AI namely, the Artificial Neural Network (ANN), Fuzzy 
Logic (FL) and Evolutionary Computation (EC) [8]. The three main branches of AI have their own different 
applications.On the other hand, EC is commonly used for population-based metaheuristic optimization 
process. Generally, there are four main branches for Evolutionary Computation (EC) which are Genetic 
Algorithm (GA), Genetic Programming (GP), Evolutionary Strategy (ES) and Evolutionary Programming 
(EP) [9]. In addition, EP and DEA are categorized as the population-based techniques. These meta-heuristic 
methods are also classified as part of CL 

This paper presents a Hybrid Standalone Photovoltaic System sizing optimization based on load 
profile. There are three optimization techniques have been developed, used to determine the optimum PV 
modules, charge controllers, inverters, battery and diesel generator. Comparative studies have been 
performed PEP, CEP and DEA; while, the ISA used as the benchmark technique. 


2. SYSTEM DESCRIPTION 

This section illustrates the method for sizing HSAPV system. The HSAPV system seen in this 
research is a PV-battery-ganset power system with AC load as shown in Pigure 1. 



Pigure 1. Block diagram of HSAPV system 


The main thing before forming a HSAPV system, load demand should be determined first. In this 
research, data were collected from a rural area in Kalabakan, Sabah, Malaysia. The estimated daily load 
demand and monthly solar radiations are shown in Table 1 and Table 2 respectively. 


Table 1. Daily Load Profile 


Appliance 

Number of 
units 

Power 
per unit 

Power 

factor 

Usage 

Times 

Energy 

Usage 

Times 

Energy 

Surge 

factor 






Dry season 


Wet season 







(jan..jun) 


(juL.dec) 


240 ac, 50Hz 


w 

w 

h 

Wh 

h 

Wh 


Water pump 

3 

120 

0.7 

2 

720 

1 

360 

5 

Incandescent 

bulb 

15 

40 

1 

10 

6000 

13 

7800 

1 

Television 

3 

80 

0.95 

8 

1920 

6 

1440 

1 

Stand fan 

3 

60 

0.6 

4 

720 

10 

1800 

3 


Daily Average AC Load Demand (Wh) 


9360 


11400 



BEEI, Vol. 7, No. 2, June 2018 : 153 - 160 

























BEEI 


ISSN: 2302-9285 


□ 155 


Table 2, Expected Mothly Solar Irradiation [11] 


month 

jan 

feb 

mar 

apr 

may 

jun 

jul 

aug 

sept 

oct 

nov dec 

Irradiation, 

G (Wh/m") 

4000 

4320 

4400 

4460 

4500 

4700 

5100 

5180 

5100 

5090 

5110 5050 


2.1. System Determination 

The steps for sizing the SAPV system are defined as follows [10]: 

Step 1: Determine the required total energy demand, Erequired_dmiy 

F — rn 

^required_daily 

"inv 

Where Eac is the total energy required daily for AC load as delivered by the battery in Wh, 77 , 01 /is the average 
efficiency of the inverter used for the designed AC loads. 

Step 2: Calculate the System Voltage (SV) for the HSAPV system. The SV can be selected as 12V, 24V or 
48 V depending on size of load demand required as shown in Table 3. 


Table 3. Range ofSV[10] 



Condition 

Recommended 

SV (V) 

Case 1 

Erequired_daily^ 1 kWh 

12 

Case 2 

1kWh < Erequired_daily< 4 kWh 

24 

Case 3 

Erequired_daily — 4 kWh 

48 


Step 3: Determine the total capacity required daily, Crequiredjiaiiy 


^ _ ^required_daily 

^required_daily ~ 


Step 4: Determine the battery bank capacity required, Cbank_required 


^ ^ T^autonomy 

^bankjrequired ^required^daily ^ dOD 


( 2 ) 


(3) 


where DOD^ax is the maximum battery depth of discharge. Tamonomy represents the expected number of days 
that the load demand will be happened by battery bank without sunlight availability. 

Step 5: Determine the total load current from the battery bank, Itotaijoad_current and the battery bank discharge 
rate, 


totalJoad_current 


SV 


X DCPower -I- 


X ACPower 


^inv 


‘ bank disch 


^bankjrequired 
^ total_load_current 


(4) 

(5) 


where Thank_disch is the discharge rate of battery bank in h. 

Step 6 : Determine the battery bank configuration, i.e. the number of battery in series string, Nseries_bank, 
number of battery strings in parallel, and bank capacity of the selected battery, Ctank_seiected, 


N. 


SV 


seriesjyank 


^ nomjyatt 


(6) 


where Vnom_batt is the nominal battery voltage and Cper_banery is the AH capacity which can selected from the 
datasheet 


f^j _ ^revised_bank_required 

^^parallel_bank ~ ^ 

^per_battery 


(7) 


-‘bank^selected 


~ ^parallel__^bank ^ ^perj)attery 


( 8 ) 
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Step 7: Based on MPPT charge controller; determine the maximum and minimum number of series PV 

module based on open circuit voltage. N s^yi^s_max_based_on_Voc tmd ^senes_mm_based_on_Voc. 


N. 


series max_baseci on Voc 


0.95 X V'l^ax.cc 


V'r 


max_oc 


(9) 


N. 


series min_based on Voc 


1.1 X V^niiii_wm_cc 
^min.mp ^ fpv_cc 


( 10 ) 


where V^ax_cc is the maximum input voltage rating of the charge controller and V^ax_oc is the maximum open 
circuit voltage rating of PV module. Vmm_wimhw_cc is the minimum allowable window input voltage to the 
MPPT of the charge controller, is the minimum voltage at maximum power of PV module and^^ ^c is 

cabling efficiency from PV to charge controller. 

Step 8: Determine the maximum number of series PV module based on maximum power voltage, 

^max_series_hQsed_on_Vmp 


N. 


0.95 XV, 


max _win_cc 


series_max_based_on_Vmp 


( 11 ) 


Where Vmax_wm_cc is the maximum window voltage of the charge controller and V„ax_mp is the maximum 
voltage at maximum power at PV module. 

Step 9: Determine the corrected output power of module, Pmod_corrected sub-system efficiency, Pmod_corrected 


Pmod_corrected Pmpjtc ^ ftemp ^ fmm ^ fdirt ( 12 ) 

Where P^p^stc is the rated power at STC, is the derating factor for dirt, is the derating factor of the 
manufacturing tolerance and/,g„,pis the derating factor for temperature. 


1; 


pv_ss 


1; 


pv_ccible 




controller 


Xii. 


batt 


(13) 


Where t]py_i;at)ie is the efficiency of cable from PV array to battery, tjcontroiier is the efficiency of the MPPT 
charge controller and is the efficiency of the battery. 

Step 10: Determine the total number of PV modules for the PV array, Njjnppt 


N 


Tjnppt ~ 


^required_daily ^ fo 
^mod_corrected ^ PSH X 


(14) 


Where 77pi,_ss is the sub-system efficiency of PV array and Pmp_stc is the rated maximum power at Standard 
Test Conditions (STC). 

Step 11: Determine the number of parallel PV modules, Np_mpp,, 


N, 


Pjnppt 


NTjnppt 

Nsjnppt 


(15) 


where Ns_mppt is the total number of PV modules in series between Nserksjnax and Nseries_mm 
Step 12: Determine the peak capacity of array using the MPPT charge controller, Parmy_mppt_cc 

Pdrrdy_mppt_cc Pmpjnod^stc ^ ^sjnppt ^ ^Pjnppt ( 

Step 13: Determine the apparent power of the inverter required for 30 minutes at maximum demand. 


^invJOmin ~ ^mdX^AC_demdnd ^ ^finv 

Step 14: Determine the apparent power of the inverter required for surge demand. 

^invjurge ~ ^mdx_AC_surge ^ ^finv 

d' _ Preal yn 

'^m.ax_AC_surge ~ PF 


(17) 


(18) 

(19) 
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where Sinvjomin is the apparent power of the inverter at maximum demand, Sinv_surge is the surge rating of the 
inverter, S„ax_Ac_surge is AC the surge load demand, S„ax_Ac_demand^s, the maximum AC load demand, Sf„yis the 
safety factor of inverter, is the real power consumption when running, Sf is the surge factor and PF is the 
power factor. 

Step 15: Determine the minimum size of diesel generator; 

c _ i^max_AC_demand~^Lnv_30mLn) /'nrw 

^genjmax^demand ~ ~ 


where Sg,.„_max_demandj,araUei is the minimum diesel generator rating based on power demand, S,„ax_Ac_demand is 
the maximum AC load, S,„v_i 0 mm is the 30 minutes’ apparent power rating of the inverter, fgg is diesel 
generator factor and/derate is the total derating factor. Sgen_max_siirge can be calculated using Equation (21). 

_ i^max_AC_surge~^inv_surge) \ 

^gen_max__surge ~ ^Jgo 


where Sge„_max_surge is the minimum diesel generator rating based on surge condition, Smax_Ac_surge is the AC 
surge load demand and Si„y_surgei& the surge rating of the inverter. 


P _ i^batt_charger ^max_AC_demand) 

^aen^batt - jyyyy 




go 


( 22 ) 


where Sgenj,att is the diesel generator rating based on battery charging and maximum load profile. 

Step 16: Calculate the Performance Ratio, PR 

PR =-- (23) 

aarray_stc ^ 

where Esys_exp is the expected kWh output of the system and Parray_stc is the nominal PV array power at STC. 
PSH represents the estimated annual PSH derived from Table 2. 


3. DOLPHIN ECHOLOCATION ALGORITHM 

Natural resources have inspired many of the man-made technologies. Researchers believe that 
dolphins are the second to humans in smartness and intelligence. Dolphin can use the echolocation method to 
help them navigate better in underwater such as to find foods or communicating. The sounds which produced 
from a dolphin are in the form of clicks which has higher frequency than the sound for communication. 
These clicks create sound waves that able to travel quickly and when the sound waves bounce off objects 
around them, they return to the dolphins as echoes and the dolphins can determine the distance [12]. 

Basically, the DEA-based optimization process, it starts with the initialization of an objective 
function. After that, the Predefined Probability of a loop, PP (LoopJ and the Accumulative Eitness of the /,/, 
alternative to be chosen for the j,h variable, AEy are calculated using 

PPiLoopd = PPdoopA + (1 - PPdoopA) X LyZZ-J <“> 

AFij = ~ X (^e ~ 1^1) ^ Fitness (t) -I- AFij_y (25) 

where PP(Loopi) is the convergence factor of the first loop. Loop, is the number of current loop and 
LoopNumber is the number of loops. AFij.j is the previous AE of the i,h alternative to be chosen for the jth 
variable. Besides, Rg is the valuable radius where the AE of the i,h alternative to be chosen for the jth variable 
is affected from its fitness. This radius is recommended to be not more than quarter of the search space [13]. 
k is the positive and negative value for the effective radius where the AE of the ith alternative to be chosen for 
the jth variable are affected from its fitness and Fitness (i) is the fitness of the location i. 

Next, a small value of s is added to AE are consistently in the search space for the selected option. 
Up to now, previous studies have indicated that, the value of e is important and superior to be less than the 
minimum value achieved for the fitness [13]. 

AF = AF + E (26) 

After that, set the AF=Q for the best candidates from the population while the rest of the candidates 
are maintained with the original value for their AF. The step can be defined as follows: 

Eor j=l: Number of variables ; i=l: Number of alternatives; If /=The best location (j) 
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AFij = 0 


Next step is to determine the probability, Py of choosing the i,h alternative for all the number of 
variables, j is calculated using: 


After that, PP(Loopi) is allocated to all alternatives chosen for all variables. Then, the next location 
of candidate is determined according to the probabilities assigned to each alternative. In this study, the PR 
value was determined based on the fitness value of each candidate. The results were ranked accordingly to 
their fitness values with the candidate having maximum PR. Finally, all steps are repeated until the maximum 
loop numbers is reached. 



Figure 2. Flowchart of Sizing and optimization technique using DBA 
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4. RESULTS AND ANALYSIS 

This section describes the results of HSAPV system using DEA and EP models with ISA as shown 
in Table 4. These results finding the maximized PR and the optimal configuration among set of system 
components which can meet the load demand requirements. Prom the table, it is observed that all the four 
technique ISA, DEA, CEP and PEP conclude no of inverter code is 7. Beside that, the number of charge 
controller code identified by ISA, DEA and PEP is identical, i.e. 7; while, CEP determines 8 as the charge 
controller code value. ISA is the technique for the benchmarking of the whole optimization algorithm 
techniques. Apart from that, the performance DEA and PEP was found to produce the optimal value of PR 
with ISA, i.e. 0.7682. However, DEA has speed up the computation time to accomplish the minimal value. 
This result shows 91.32% faster than ISA technique. 


T ab le 4. Performance of Comparison of DEA and EP Models with I S A 


Results 

ISA 

EP Techniques 
CEP FEP 

DEA 

Battery code 

12 

12 

12 

12 

PV code 

12 

6 

12 

12 

Inverter code 

7 

7 

7 

7 

Charge Controller code 

7 

8 

7 

7 

Diesel Generator Code 

11 

11 

11 

11 

N_batt_s 

8 

8 

8 

8 

N_batt_p 

2 

2 

2 

2 

N_pv_s 

2 

3 

2 

2 

N_pv_p 

8 

6 

8 

8 

Nt inv 

1 

1 

1 

1 

Nt_cc 

2 

2 

2 

2 

Nt_diesel_gen 

1 

1 

1 

1 

Optimal, PR 

0.7682 

0.7428 

0.7682 

0.7682 

Comp, time, (s) 

2003.04 

174.58 

174.02 

173.91 


The next analysis is the maximization of PR based on four techniques by using different population 
size is tabulated in Pigure 3. Based on figure, for the small size population (n=5), the range of maximum PR 
for three optimization is 0.71-0.750 compare with ISA (i.e 0.7628). Purthermore, when the population size is 
increased, the values of fitness values are increase accordingly. On the other hand, the PEP and DEA reached 
the maximum value of PR (i.e 0.7628) with population size of 30 as those set by ISA, while the CEP is 
reached at value 0.7448 until the end of population size. In addition, the DEA is maintained with the 
maximum value of PR at population size of 25 onwards (i.e 0.7628), while the PEP is reached the maximum 
PR at population size of 30. These mean that DEA is better than CEP and PEP in optimizing the HSAPV 
system for maximization of fitness value. 



» ISA 
■ CEP 
A FEP 
DEA 


Pigure 3. Performance of fitness value with different techniques at different number of population 


5. CONCLUSION 

This paper has presented Hybrid Standalone Photovoltaic Systems Sizing Optimization based on 
load profile.There are three different optimization techniques were investigated before selecting the best 
optimal value based on sizing algorithm. Prom the results, the DEA has reached the maximum of fitness 
values, together with the minimal computation time. The DEA has been compared with the ISA, whereas 
ISA is the iterative technique where can produced the optimal value of PR. The development of DEA is also 
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feasible for solving the other optimization problems. For future works, DBA could be potentially conducted 
to apply in multisource of SAPV systems. 
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